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Abstract. The Trojan Horse Method was applied for the first time to a Radioactive Ion Beam 
induced reaction to study the reaction ^**F(p,q)^^0 via the three body reaction ^**F(d,Q ^^0)n 
at the low energies relevant for astrophysics. 

The abundance of ^*F in Nova explosions is an important issue for the understanding of this 
astrophysical phenomenon. For this reason it is necessary to study the nuclear reactions that 
produce or destroy ^*F in Novae. ^*F(p,a)^^0 is one of the main ^*F destruction channels. 
Preliminary results are presented in this paper. 



1. Introduction 

The 7-ray emission following the Nova explosion is dominated by the 511 keV energy line, coming 
from the annihilation of positrons produced by the decay of radioactive nuclei. Among them, 
the is especially important because of its expected abundance in the Nova environment and 
because of its lifetime, that well matches the timescale for Nova ejecta to become transparent 
to 7-ray emission. To understand the Nova explosion phenomena, it is then important to 
know the rate of nuclear reactions producing and destroying ^^F. At relevant temperatures, 
the ^^F(p,a)^^0 reaction is expected to dominate by roughly a factor of 1000 [1] and it is 
the uncertainty in this reaction rate that gives the main nuclear contribution to the overall 
uncertainty in the final abundance of ^^F. The cross section of the reaction is dominated 
by the presence of several levels in the ^^Ne compound nucleus around the ^^Fp threshold. 
Moreover, the tails of some subthreshold levels are expected to contribute to the resonance cross 




Figure 1. The THM basic diagram. The upper vertex describes the virtual decays of particle 
a into X s, while the lower one refers to the process of interest, namely b + x ^ c + C. In this 
paper, a is deuteron going into p and n and the reaction of interest is ^^F-\-p — t-^^O+q. 



section. Especially important is not only the width of these resonance states but also their 
interference effects. It is then important to explore the behavior of the cross section in the 
center of mass energy range from up to 500 keV to get information about the reaction rate at 
Nova temperatures. 

In spite of several and long lasting attempts to measure the cross section of ^^F(p,a)^^0 
reaction (see e.g. [SHUJ) the situation is still not satisfactory. In order to get new complementary 
pieces of information on this process, we performed two experiments aiming at the measurement 
of the ^^F(p,a)^^0 using two different techniques: the direct thick target method and the indirect 
Trojan Horse Method (THM). We report here only on the latter one. The THM was applied 
to a reaction induced by a Radioactive Ion Beam (RIB) for the first time in the experiment 
described here. Namely, the three body reaction ^^F(d,a ^^0)n was used to infer information 
on the process of interest, ^^F(p,Q;)^^0, in the low energy region important for astrophysics. 

2. The Trojan Horse Method 

The Trojan Horse Method (THM) [181122] has been proposed and exploited to study the cross 
sections of nuclear reactions between charged particles at astrophysical energies. Briefly, a re- 
action of the type 

x + b^c + C (1) 
can be studied by using a three-body reaction of the type 

a + b^ c + C + s, (2) 

where a is a nucleus with a strong cluster structure of the type x ^ s and it is often called the 
trojan horse nucleus. The measurement of the reaction (2) is performed at energy higher than 
both the Coulomb barrier and the energy required to breakup particle a into its components x 
and s. 

If reaction (2) proceeds through a quasi-free reaction mechanism, then it can be described by 
the diagram shown in figure [T] and the cross section of the three-body reaction can be factorized 
in three parts: a kinematical term, a term describing the virtual decay of a into x and s and 
one giving the cross section of the reaction (1). Using a Plane Wave Impulse Approximation 



approach to describe the diagram in figure [U this can be written as |19ft21j 



"'"^ oc(A-F)|*WPf^f (3) 



where KF is the kinematical factor, describes the impulse distribution of the spectator 

s inside a and ps is the momentum of the spectator. The symbol N for the cross section of 
the reaction (1) indicates that it is a pure nuclear cross section, whitout effects coming from 
Coulomb interactions. 
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Figure 2. Scheme of the experimental set-up used in the experiment. The PPAC and MCP 
detectors were used for beam tracking, hence allowing for the reconstruction of the incidence 
angle of the beam on the CD2 target (the oval in the figure) The reaction products are 
then detected by the plane made of bidimensional position sensitive detectors and by the two 
multistrip detectors. The safety disk protects the strips at very low angles from the risks that 
can derive from being hit by the direct beam. 

Hence, by measuring the cross section of the three body process (2) and the impulse 
distribution of s inside a (represented by |<^(ps)p), and by calculating the kinematical factor 
(KF), one can derive the cross section for the process of astrophysical interest (1). The cm. 
interaction energy of the reaction (1) is given by the relation Efjx = Ecc — Q2, where Ecc is the 
relative energy of the outgoing particles c and C and Q2 is the q-value of the two body process 
(1). The calculated two body cross section is a purely nuclear one, the Coulomb barrier having 
been overcome by the TH mechanism |18H2m[^^l26j . Hence, in order to compare the THM cross 
sections with the directly measured ones, the barrier suppression effects have to be taken into 
account by correcting the THM cross section for the probability of penetration under barrier. 

The THM has been used a number of times [TH1E3H33] to measure cross sections of 
astrophysical interest and it is presently regarded as a very powerful tool by the nuclear 
astrophysics community. Moreover the THM can be also used to study neutron induced 
reaction [M] . 

The simple theoretical approach discussed above has been improved along the years. More 
sophisticated techniques have been introduced to treat more complicated situations as in the 
cases where the cross section is dominated by the presence of narrow, close and/or interfering 
resonances [28l[29]. 



3. The Experiment 

The beam was produced at the CRIB separator of the Center for Nuclear Study (CNS) of 
the University of Tokyo installed at RIKEN campus in Wako, Japan, by using the ^^0(p,n)^^F 
reaction. The ^^0 primary beam (8^, E=4.5-5 MeV/A) was provided by a AVF cyclotron. It 
impinged on a primary gas target of H2. The reaction products were then separated by the 
CRIB doubly achromatic spectrometer, described in detail in ref. [35]. 
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Figure 3. Beam image on the PPAC detector (a), on the MCP detector (b) and beam spot on 
target position reconstructed event by event (c). 

Two beam production tests were performed before the experiment. The characteristics of 
the ^*F beam used during the experiment were the following: beam energy peaked at 47.9 MeV 
(FWHM 1.9 MeV) with maximum beam intensity of 2x10^ pps and purity better than 98%. 
The intensity was higher than 5x10^ pps throughout the measurement. This beam was used to 
bombard a CD2 secondary target. 

The experimental set-up, schematically shown in figure was very different from the typical 
ones used for THM experiments. Indeed, it consisted essentially of two parts: a beam tracker 
placed before the secondary target and two planes of position sensitive silicon detectors to detect 
the ^^O and a ejectiles. 

The tracks of the beam particles were reconstructed, event by event, using two Parallel Plate 
Avalanche Chambers (PPAC) with a position resolution (x-y) of roughly 1 nini and a distance 
between the two planes of 500 mm. During the experiment the second PPAC broke down and 
it was replaced by a MicroChannel Plate detector (MCP) having the same position resolution. 
The MCP detected the electrons coming from the interaction of the beam with a Carbon foil 
placed in the same position of the broken PPAC. The precise particle tracking allowed for the 
determination of the intersection point of the projectiles on the secondary target. This is a key 
point for the application of the THM, as discussed in [36]. 

The secondary target was made of a thin (roughly 150 /ig/cm^) CD2 foil of 30 mm diameter. 

The detector setup was based on the modular system built at the INFN Laboratori Nazionali 
del Sud called ASTRHO (A Silicon Array for TRojan HOrse). In the configuration used for this 
experiment, ASTRHO consists of a mechanical frame that hosts two planes of position sensitive 
detectors. The first plane consisted of 8 bidimentional position sensitive silicon detectors 
manufactured by Hamamatsu. Each detector is 45x45 mm^ and measures the energy and 
the position on x and y directions of the impinging particle. The spatial resolution is of 1 mm in 
both directions. The second plane was made of a pair of double sided silicon multistrip detectors 
(DSSSD), manufactured by Micron. Each detector (59x50 mm^) has 16 xl6 strips in x and y 
direction respectively, with a spatial resolution of each strip of 1.5 mm. 



The detector holders were especially designed for this new apparatus, with the aim of 
having a well defined geometry by construction and a support that can be adapted for different 
experiments by changing the distances among the different detector planes and the target. 

In the present case the first plane of bidimensional position sensitive detectors was optimized 
to detect the a particles, while ^^O is detected at most forward angles by DSSSD. The angular 
ranges covered by the experimental set-up were 11° < 9 < 31° for a particles and 2° < 9 < 12° 
for ^^O ions. For each detector the x and y position and the energy of the impinging particles 
were recorded together with the time of flight. 

4. Data Analysis 

In order to apply the THM to the study of ^^F(p,a)-^^0, we selected the ^^F(d,a ^^0)n three 
body reaction. 

Given the low intensity of a RIB with respect to a stable beam, the use of a large solid 
angle setup like ASTRHO was necessary to increase the global detection efficiency. Moreover 
the beam track reconstruction event by event was mandatory to reconstruct the angles of the 
ejected particles with respect to the actual incident particle trajectories, improving the angular 
resolution of the experiment. 

In figure E] the beam image reconstructed on the two detectors planes and on the target are 
reported. It is clearly seen that even if the beam spot on the target position is small, the beam 
divergence is not negligible. 




Figure 4. The Q-value spectra for the events that pass the conditions imposed to identify the 
three body reaction channel of interest. The arrow represents the position for the theoretical 
Q-value of the ^^F(p,a)^^0 reaction (0.658 MeV). The hue is the result of a gaussian fit and the 
obtained parameters are reported in the figure. 

The calibration in energy of the detectors was performed by using the elastic scattering peaks 
measured at different beam energies for the reaction ^®0-|-^^^Au. Moreover, a triple alpha source 
was used to calibrate at low energies. The detectors made by Hamamatsu were calibrated in 



position by using grids placed in front of each detector during the cahbration runs. The energy 
lost on the target by the beam, the detected particle energy loss on the target and on the dead 
layer of the detectors were also taken into account for the calibration in position and energy. 

After calibration, the data analysis of a THM experiment mainly proceeds through the 
following steps: identification of the process of interest, selection of the events corresponding to 
the quasi- free reaction channel, extraction of the two body cross section of astrophysical interest. 

In the present experiment, due to the limitation in mass and Z identification of the outgoing 
particles, other three reactions with three bodies in the final state can give contribution in 
the phase space region where one expects events coming from ^^F(d,a ^^0)n. These reactions 
are ^^F(d,a ^^N)p, ^SF(d,p ^^0)p and ^^F(d,p ^^F)n. We were able to select the channel of 
interest by applying convenient cuts in the phase space. Figure 4 shows the Q-value spectra 
that was obtained by applying the event selection procedure. The position and width of the peak 
are in very good agreement with the experimental conditions. The expected value for events 
coming from ^^F(p,a)^^0 is 0.658 MeV to be compared with a measured one of 0.647 MeV. 
Also, the measured width of the Q-value peak of FWHM 2.3 MeV is well accounted for by the 
intrinsic beam energy spread (1.9 MeV FWHM) plus the kinematical enlargement. This gives 
us confidence of the correct identification of the events coming from the channel of interest. 
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Figure 5. Momentum distribution for the p-n intercluster motion in deuteron. The solid line 
is the Hulthen function in momentum space. 

The second step mentioned above in a THM analysis is the selection of events that come 
from the quasi-free channel in the ^^F(d,a ^^0)n reaction, among other processes falling the 
same particles in the final state. The strongest evidence of the predominance of the quasi- free 
mechanism is given by the shape of the momentum distribution for the p-n intercluster motion 
in deuteron. This is shown in figure by black dots as a function of the momentum of the 
spectator particle, ps- 

The solid line in figure [S] represents the Hulthen function in momentum space with the 
standard parameters values [37]. This is expected for deuteron as trojan horse nucleus if the 
plane- wave impulse approximation is used to describe the reaction mechanism. Distortion effects 



influence the behavior of the distribution at high values of momentum, so the rest of the analysis 
was done by imposing a cut on the momentum of the spectator particle. Namely, only events 
where the spectator momentum was lower than 50 MeV/c were accepted. 
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Figure 6. The nuclear cross section spectrum in function of the p-^^F cm energy for the events 
that pass the conditions described in the text. 

Assuming the events selected according to the previous procedure are actually coming from 
the quasi-free contribution to the reaction yield, then one can apply equation (3) to obtain the 
cross section of interest for the ^^F(p,Q!)^^0 process at astrophysical energies. The result is 
shown in figure 6, that represents the nuclear cross section for ^'^F(p,a)^^0 measured by THM 
down to zero (and even negative) energies without any suppression due to the Coulomb barrier. 
In order to compare the present result with the direct data, it is necessary to correct the obtained 
nuclear cross section for the probability of penetration of the Coulomb and, possibly, centrifugal 
barrier. Data analysis is still in progress. 

5. Conclusion 

The THM was applied for the first time to study a reaction induced by a radioactive ion beam. 
Even with the use of indirect methods like THM, the measurement of cross sections of interest 
for nuclear astrophysics remains one of the most difficult tasks in nuclear physics, especially 
if one has to use RIBs when the low beam intensity add on the top of the low cross sections 
typical of astrophysical nuclear processes. The use of a new large solid angle detector setup 
especially designed for this experiment was hence necessary in order to recover, at least in part, 
the possibility of accumulating a statistics high enough to study the processes of interest. 

The preliminary data analysis has demonstrated the feasibility of the experiment and gives 
us confidence that the results coming from the final analysis of the data will be helpful in 
understanding the problem of ^^F destruction in Novae. 

6. Acknowledgments 

We are indebted to the AVF and CRIB technical staff for their invaluable contribution to the ex- 
periment. Also we would like to thank the INFN-LNS workshop crew headed by Mr.B. Trovato 




for the perfect construction of the ASTRHO mechanical frame and Mr. C. Marchetta from 
the target laboratory at INFN-LNS for providing us with, among others, high quality and high 
purity CD2 targets. The work was supported in part by the Italian Ministry of University and 
Research under Grant No. RBFR082838 (FIRB2008). 

[1] A. Coc et al., A.&A., 357 (2000) 561 and references therein. 

[2] M. Wiescher and K.U. Kcttner, Astrophys. J. 263 (1982) 891. 

[3] K.E. Rehm, M. Paul, A.D. Roberts, et al., Phys. Rev. C52 (1995) R460. 

[4] K.E. Rehm, M. Paul, A.D. Roberts, et al., Phys. Rev. C53 (1996) 1950. 

[5] K.E. Rehm, C.L. Jiang, M. Paul, ot al., Phys. Rev. C55 (1997) R566. 

[6] J.S. Graulich, F. Binon, W. Bradfield-Smith, et al., Nucl. Phys. A626 (1997) 751. 

[7] R. Coszach, M. Cogncau, C.R. Bain, et al., Phys. Lett. B353 (1998) 184. 

[8] S. Utku, J.G. Ross, N.P.T. Bateman, et al., Phys. Rev. 057 (1998) 2731 and C58 (1998) 1354. 

[9] Y.M. Butt, J.W. Hammer, M. Jaeger, ct al., Phys. Rev. C58 (1998) RIO. 

[10] J.-S. Grauhch, S. Chcrubini, R. Coszach, ot al.,Phys. Rev. C63 (2001) 011302 and references therein. 

[11] D.W. Bardayan, J.C. Blackmon, W. Bradfleld-Smith, et al., Phys. Rev. C63 (2001) 065802. 

[12] D.W. Bardayan, J.C. Batchelder, J.C. Blackmon, et al., Phys. Rev. Lett. 89 (2002) 262501. 

[13] N. De Sereville, A. Coc, C. Angulo, et al., Phys. Rev. C67 (2003) 052801 (R) and references therein. 

[14] R.L. Kozub, D.W. Bardayan, J.C. Batchelder, et al., Phys. Rev. C71 (2005) 032801 (R). 

[15] R.L. Kozub, D.W. Bardayan, J.C. Batchelder, ct al., Nucl. Phys. A 758 (2005) 753c-756c. 

[16] C. E. Beer et al., Phys. Rev. C 83, 042801 (2011) 

[17] A. S. Adekola, D. W. Bardayan, J.C. Blackmon, et al., Phys. Rev. C83, 052801 (R), (2011). 

[18] C. Spitalcri. in Proceedings of the Fifth Hadronic Physics Winter Seminar, Folgaria - TN, Italy, Ed. World 
Scientific, Singapore, 1990. 

[19] S. Chcrubini, V. N. Kondratyev, M. Lattuada, C. Spitaleri, D. Miljanic, M. Zadro, and G. Baur. ApJ, 457, 
855, (1996) 

[20] G. Baur and S. Typel. Prog. Theor. Phys. SuppL, 154, 333, (2004) 

[21] S. Typcl and G. Baur. Ann. of Phys, 305, 228, (2003) 

[22] A. M. Mukhamedzhanov et al. Nucl. Phys., A787, 321C, 2007. 

[23] C. Spitaleri, M. Aliotta, S. Chcrubini, M. Lattuada, Dj. Miljanic, S. Romano, N. Soic, M. Zadro, and R. A. 

Zappala, Phys. Rev. C60, 055802, 1999. 

[24] C. Spitaleri et al, Phys. Rev. C63, 005801, 2001. 

[25] C. Spitalcri et al, Phys. Rev. C69, 055806, 2004. 

[26] A. Tumino et al., Phys. Rev. Lett, 98, 252502, 2007. 

[27] M. La Cognata et al., Phys. Rev. C76, 06584, 2007. 

[28] M. La Cognata et al., Phys. Rev. Lett, 101, 152501, 2008. 

[29] M. La Cognata et al., ApJ, 708, 796-811, 2010. 

[30] Adelberger, E.G., et al. Solar fusion cross sections. II. the pp chain and CNO cycles (2011) Rev. Mod. Phys., 

83 (1), pp. 195-245. 

[31] M.L.Sergi et al, Phys. Rev. C82, 032801 (2010) 

[32] L.Lamia et al,J. Phys. G: Nucl.Part.Phys. 39, 015106 (2012) 

[33] S. Cherubini et al,Act Phys. Pol 5 42,769 (2011) 

[34] M. Gulino et al,J. Phys. G: Nud.PartPhys. 37, 125105 (2010) 

[35] Y. Yanagisawa, S. Kubono, T. Teranishi, et al., Nucl. Instr. Meth. A 539 (2005) 74 

[36] S. Cherubini et al,AIP Gonf Procee. proceedings of OMEGll, in press (2012) 

[37] M. Zadro et al, Phys. Rev. C40 181, 1989 



